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Abstract 

The aim of this paper is to study direct methods for the analyticity 
of the pressure for certain classical unbounded spin models. We provide 
a representation in terms of the Witten Laplacian on one-forms of the 
nth-derivative of the pressure as function of some order parameter t. The 
technique involves the formula for the covariance introduced by B. Helffer 
and J. Sjostrand. 

1 Introduction 

As already mentioned in [66] , The methods for investigating critical phenomena 
for certain physical systems took an interesting direction when powerful and so- 
phisticated PDE techniques were introduced. The methods are generally based 
on the analysis of suitable differential operators 



and 



W^^^ = -A+^-^ —+ Hess$. 



These are in some sense deformations of the standard Laplace Beltrami oper- 
ator. They are commonly called Witten Laplacians, and were first introduced 
by Edward Witten, [18] in 1982 in the context of Morse theory for the study 
of topological invariants of compact Riemannian manifolds. In 1994, Bernard 
Helffer and Johannes Sjostrand [8] introduced two elliptic differential operators 



and 



A^°^ := - A + V$ • V 



A';^'^ := -A + V$ • V + Hess$ 
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sometimes called Helffer-Sjostrand operators serving to get direct method for 
the study of integrals and operators in high dimensions of the type that ap- 
pear in Statistical Mechanics and Euclidean Field Theory. In 1996, Johannes 
Sjostrand [13] observed that these so-called Helffer-Sjostrand operators were in 
fact equivalent to Witten's Laplacians. Since then, there have been significant 
advances in the use of these Laplacians to study the thermodynamic behavior 
of quantities related to the Gibbs measure Z~^e~'^dx. 

Numerous techniques have been developed in the study of integrals asso- 
ciated with the equilibrium Gibbs state for certain unbounded spins systems. 
One of the most striking results is an exact formula for the covariance of two 
functions in terms of the Witten Laplacian on one forms leading to sophisti- 
cated methods for estimating the correlation functions of a random field. As 
mentioned in [10], this formula is in some sense a stronger and more flexible 
version of the Brascamp-Lieb inequality [1]. The formula may be written as 
follow: 

cov(/, g)= J (4'^" V/ • Vg) e-^^^Ux. (1) 

We attempt in these notes, to study a direct method for the analyticity of 

the pressure for certain classical convex unbounded spin systems. It is central 
in Statistical Mechanics to study the differentiability or even the analyticity 
of the pressure with respect to some distinguished thermodynamic parameters 
such as temperature, chemical potential or external field. In fact the analytic 
behavior of the pressure is the classical thermodynamic indicator for the absence 
or existence of phase transition. The most famous result on the analyticity of 
the pressure is the circle theorem of Lee and Yang [28] . This theorem asserts 
the following: consider a { — 1,1} —valued spin system with ferromagnetic pair 
interaction and external field h and regard the quantity 2; = as a complex 
parameter, then all zeroes of all partition functions (with free boundary condi- 
tion), considered as functions of z lie in the complex unit circle. This theorem 
readily implies that the pressure is an analytic function of h in the region h> 
and h < 0. Heilmann [29] showed that the assumption of pair interaction is 
necessary. A transparent approach to the circle theorem was found by Asano 
[30] and developed further by Ruelle [31], [32], Slawny [33], and Gruber et al 
[34]. Griffiths [35] and Griffiths-Simon [36] foimd a method of extending the 
Lee- Yang theorem to real- valued spin systems with a particular type of a priory 
measure. Newman [37] proved the Lee- Yang theorem for every a priory mea- 
sure which satisfies this theorem in the particular case of no interaction. Dunlop 
[38] , [39] studied the zeroes of the partition functions for the plane rotor model. 
A general Lee- Yang theorem for multicomponent systems was finally proved by 
Lieb and Sokal [40]. For further references see Glimm and Jaffe [41]. 
The Lee- Yang theorem and its variants depend on the ferromagnetic character 
of the interaction. There are various other way of proving the infinite differ- 
entiability or the analyticity of the pressure for (ferromagnetic and non ferro- 
magnetic) systems at high temperatures, or at low temperatures, or at large 
external fields. Most of these take advantage of a sufficiently rapid decay of cor- 
relations and /or cluster expansion methods. Here is a small sample of relevant 
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references. Bricmont, Lebowitz and Pfistcr [42], Dobroshin [43], Dobroshin and 
Sholsman [44], [45], Duneau et al [46] , [47] , [48] , Glimm and Jaffe [41], [49], Israel 
[50], Kotecky and Preiss [51], Kunz [52], Lebowitz [53], [54], Malyshev [55], Maly- 
chev and Milnos [56] and Prakash [57]. M. Kac and J.M. Luttingcr [58] obtained 
a formula for the pressure in terms of irreducible distribution functions. 

In this present study, we propose a new way of analyzing the analyticity of 
the pressure for certain unbounded models through a representation by means 
of the Witten Laplacians of the remainder of the Taylor series expansion. The 
methods known up to now rely on complicated indirect arguments. 

2 Towards the analyticity of the Pressure 

Let A be a finite domain in Z"^ {d > 1) and consider the Hamiltonian of the 
phase space given by, 

$(a;) = $A(a;) = y + -^{x), x G M^. (2) 

where 

]a"V«'|<C„, VaeNl^l, (3) 

Hess$(a;) > So, < So < 1. (4) 

Let g is a smooth function on with lattice support Sg = T. We identified 
with g defined on by 

g{x) = g{xr) where x = {xi)-^^ and a;r = {xi)^^^' (5) 
and satisfying 

\d"Vg\ <Ccc Va e N'^I (6) 

Under the additional assumptions that ^ is compactly supported in and g is 
compactly supported in R'", it was proved in [66] (sec also [8]) that the equation 

has a unique smooth solution satisfying V'^/(a;) — > as |a;| — > oo for every 
k>l. 

Recall also that V/ is a solution of the system 

(- A + V$ ■ V) V/ + Hess$V/ = Vg in R^. (7) 

As in [66] and [8], these assumptions will be relaxed later on. 
Let 

'fi{x)=^x)-tg{x), (8) 
where x = {xi)i^A, and assume additionally that g satisfies 

Hessg < C. (9) 
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We consider the following perturbation 



eA{t) = log 



:/ 



(10) 



Denote by 



and 



/ 



dxe-*A(^) (11) 



< ■ >t,A= ^ ^ . (12) 



3 Parameter Dependency of the Solution 

From the assumptions made on <1> and g, it is easy to see that there exists T > 
such that or every t G [0, T), satisfies all the assimiptions required for the 

solvability, regularity and asymptotic behavior of the solution f{t) associated 
with the potential $:y(.t). Thus, each t G [0,T) is associated with a unique 
C°°— solution, f{t) of the equation 

Hence, 

A'i}j{t)=Vg (13) 
where v{t) = V/(f). Notice that the map 

t ^ v(t) 

is well defined and 

{v(t):ie[o,r)} 

is a family of smooth solutions on satisfying 

d"v{t)^0 as|a;|^oo Va e N'^I and for each i e [0, T) 
and corresponding to the family of potential 

{$i:te[0,T)}. (14) 



Let us now verify that v is a smooth function of f G (0,T).We need to prove 
that for each t e (0,T), the limit 

^.^v(i±£)^ 

exists. Let 
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We use a technique based on regularity estimates to get a uniform control of 
v'^(t) with respect to e. 

With £ small enough, we have 



= -A 

+ 



w{t + e)- v{t) 



£ 

Hess$*+^v(i + £)- Hess$*v(t) 



^ V$*+^ . Vvft + g) - • Vv(^) 



Equivalently, 

"v(i + e)- v(i) 



-A 



e 

■v{t + £) - ^^{t) 

£ 



^ V$*+^ • V [wit + g) - vft)] 



v(t) - : • Vv(t) 



and 



+Hess$*+' 

-( 

-Av''(f) + V$*+^ • Vv''(t) + Hess$*+^v'^(i) 



= -( )vW-( j-VvW. 

Let w(t) be the unique C°°— solution of the system 

- Aw(i) + V$* • Vw(i) + Hess$*w(t) = Hesss'v(t) - Vg ■ Vv(i). (15) 

Combining the last two systems above, we get 

-A [w(t) - v^(i)] + V$* • V [w(t) - v^(t)] + Hess$* [w(t) - v^(t)] 

'Hess$*+^ -Hess3>*~ 



= Hesss'v(t) - Vff ■ Vv(<) + v(t) 



e 



+ ( — — 1 • Vv(t) + (V$*+^ - V$*) • Vv=(t) 

+ (Hess$*+'^ - Hess$*) w^{t). 
Now using the unitary transformation U^t, we get 



(16) 



f-A+I^l^ ^j(w(i)-v^(i))e-*V2 

+Hess$* (w(t) - v^(i)) e-*'/2 (17) 
= Oe(l)e-**/2 + [(v$*+^ - V$*) ■ Vv"(i) 
+ (Hess$*+^ - Hess$*) v^(i)]e-*V2 

Next, we propose to estimate the last two terms of the right hand side of this 
equation. 
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Again using the unitary transformation U^t+s we reduce the system 



-Av^(t) + V$*+'^ • Vv^it) + Hess$*+^v^(t) 
/Hess$*+= -Hess$*' 



v(t) 



(18) 



into 



-A+ 



4 2 
/Hess$*+^ - Hess$* 



Vv(0 

+ Hess$*+'^V'^ 
v(t)e-**+V2 



(19) 



V<I>' 



Vv(t)e- 



whoro = v^(t)e *'"''^/2 Taking scalar product with on both sides of this 
last equality and integrating, we get 



t+s 



-/ 



1,2 



+ / Hess$*+^V= • V'dx 



Hess$*+^ - Hess$' 



v{t)e 



/(:^^!^!^).Vv(t)e— /2 



V^dx 



(20) 



Now using the uniform strict convexity on the left hand side and Cauchy- 
Schwartz on the right hand side, we obtain 



II II go < C for small enough e. 



We then deduce that 



(21) 



(22) 



where 



'Hess$*+^-Hess$* 



v(t)e- 



+ — - — - Hess$*+^V^ 



• Vv(t)e- 



(23) 

is bounded in uniformly with respect to s for s small enough. 

Now taking scalar product with on both sides of (28) and integrating by 

parts, we obtain 



t+s\ 



< U 



(24) 
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It follows that is uniformly bounded with respect to e in -8^*+= ^o^^ ^ small 
enough. 

Next, observe that 

||_A+r^ j V- = (25) 

where 

qe = qe - ^- ^ ' V- + ^— 1 ^ (26) 

is uniformly bounded in with respect to e for small enough e. Using regu- 
larity, it follows that for small enough £, is uniformly bounded in B\t with 
respect to e.This implies that is uniformly bounded in B\t for e small enough. 
Again, we can continue by a bootstrap argument to consequently get that for e 
small enough, is uniformly bounded in B^t for any k. 
It is then clear that for small enough e, the right hand sides of (23) is 0{e) in B^ 
and consequently, using the same argument as above, wc get that (w(f) — v'^(t)) e~* 
is 0(e) in B^t] again iterating the regularity argument, we obtain that for small 
enough e, (w(t) — v^(t)) e~*'/^ is 0(e) B^t for every k. We have proved: 

Proposition 1 Under the above on $ and g, there exists T > so that for each 
t G (0,T), v=(t) converges to w{t) in C°°. 

Remark 2 The proposition establishes that v(t) is differentiable in t and ^^(^) 

is given by the unique — solution w(t) of the system 

- Aw{t) + V$* • Vw(i) + Hess$*w(t) = Hesssfv(t) - Vg ■ Vv(i). (27) 
Iterating this argument, we easily get that, v{t) is smooth in t & {0,T). 
Now we are ready for the following: 



4 Formula for e^''\t) 

For an arbitrary suitable function f{t) = f{t, iv) 
d 



Hence, 



< fit) >t,A=< / '(i) >t,A +cov(/, g). (28) 



^ < fit) >t,A=< f 'it) >t,A + < (V/) • Vg >t,A . (29) 



Let 



A,/ := (V/) • V5. (30) 
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Thus, 

^ < /(i) >t,A=< ( ^ + ^9 ) / >*,A 



d_ 
di 



The hnear operator Ag wih be denoted by Hg. 



d 



^aW = <9 >t,A 

= < H°g >t,A; 



< ( ^+^9) 5>t.A 



d_ 
di 

d 

< 



-(4r(V5)-v5) >*,A 



By induction it is easy to see that 

g-^+Agj g>t,A (Vn>l) 

Next, we propose to find a simpler formula for (i) that only involves A 



Hgg = 4V"(Vff)-Vff 
= Agg 



Hlg = I V/ • V5 + V (V5) • V5)) ■ Vg 



where / satisfies the equation 



V/ = {Vg) . 



With v{t) = V/, as before, we get 



^ V/ • Vg = (Hess5v(<) - Vg ■ Vv{t)) ■ Vg 



and becomes 



(Hess5v(t) - Vg ■ Vv(i)) + V (^V" (Vg) • Vg)] • Vg 



4V"2V(A,5)-V<? 



Proposition 3 // 



^aW = log 



where 



is as above then 9^^\t), the nth— derivative of0A{t) is given by the formula 



and for n > 1 



9'j,{t) =< g >t,A, 



9'jC\t) = {n-l)\<A--'g> 



Proof. We have already estabUshed that 



e^jCHt)=<H^-'g>t,A forn>l. 
It then only remains to prove that 

Hr'g={n-iy.Ar'g forn>l. 

The result is already established above for n = 1,2,3,. By induction, assume 
that 



H^-'g = {n-l)\A^-'g. 



if n is replaced hy h <n. 



= (g-^+A,]{{n-iy.A--'g) 



Now 



= V(p„-Vg 



Vg 



where 

We obtain, 
d 



(i)-i / d 



A 



<i>t 



n-2. 



dt 



We then have 



d 



dt 9 



Thus, 



d_ 

dt 



( ^ + Hess5V^„ -Vg-V (V^J 



•V5 



4*"' (^V^r'5 + V(V(p„-V5) 



V5 



= A„ 



= A,H,{A--^g). 



A H 
' ' V(n-2)!"9 



^o" ^''5 ) (from the induction hypothesis) 



-A H^'^~^^a 



(n-2)! « s 

- — —TTT^g ((n — l)!^g^^5) (still by the induction hypothesis) 

(n — z). 

(n - 1)A^5. 



H^g = (n-l)!(n-l + l)A^5 



Proposition 4 If g{0) = 0, t/ien the formula 



?i")(t) = (n-l)!<A"-i3>t,A, n>2 



sii?/ /io?ds if we no longer require ^ and g to be compactly supported in R^. 
Proof. As in [8], consider the family cutoff functions 

X = Xe (34) 
(e 6 [0, 1]) in C~(]R) with value in [0, 1] such that 
X 1 for \t\ < 



\x^'Kt)\<c,^ 



for k G 
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We could take for instance 

xM = m^\t\) 

for a suitable /. 

We then introduce 

^,{x)=Xei\x\)^, arGM^ (35) 

and 

g,ix) = Xeil^Da xeR^ (36) 

One can check that both "if^ (x) and ge {x) satisfies the assumptions made above 
on and g. Now consider the equation 

- A/, + V$* • V/, = <?,- < g, >t,A. (37) 

which implies 

(- A + Vi>* • V) (8) Ve + Hess$* = V^^ (38) 

where 

It was proved in [8] that = A^^t converges in C°° to A^^t as 

£^0. ■ 

Remark 5 // we denote by Rn the remainder of the Taylor series expansion of 
the pressure PA{t), given by 

we have 

_ Pi-^'\t,) 
^" ~ (n + 1)! 

< A];g >t,A 



(n + l)|A| 

// $ and g are such that < A^g >t,K is uniformly bounded with respect to n 
and does not grow faster than |A|, we automatically get the analyticity of the 
pressure in the thermodynamic limit. 

5 Some Consequences of the Formula for nt/i— Derivative 
of the Pressure. 

In the following, we shall additionally assume that 

Vg{0) = 0, and 
V$^(0)=0 for alH e [0,T). 
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When n = 1, we recall that A'^g = g, 

^Ait) =< 9 >t,A 

and if 

v(t) = V/ = Vg, 

then we have 

(-A + V$5v • V) (g) v{t) + Hess$5^v(t) = Vg 

and as in [8] v{t) is a solution of the equation 

g=<g >t,A +y{t) ■ - divv{t). (39) 

Using the assumptions above, we have 

^aW = <9 >t,A 
= divv{t){0). 

Similarly, the formula 

= - 1)! < >t,A, 

implies that 

<^W = ("-l)!d*^v„(i)(0), 

where 

v„(i)=4Tv(A^-i5). 

Remark 6 The idea of representing 9'j^{t) in terms of diw{t){0) is due to Helf- 
fer and Sjdstrand [8] in the context of proving the exponential convergence of 
the thermodynamic limit in the one dimensional case. 

Wo conclude these notes by a discussion about the potential contribution of 
this results towards solving the two dimensional dipole gas problem. The dipole 
gas and other gases of particles interacting through Coulomb forces are very 
important statistical systems. In particular, for dipole gas, the lack of screening 
is well known [59] , and the analyticity of the pressure in the high temperature 
and low activity region has been proved in an indirect way, by means of renor- 
malization group methods (sec [60] and [61]). 

A direct proof of the analyticity of the pressure based on estimating the coeffi- 
cients of the Mayer (Taylor) series is still an open problem. The close relation- 
ship between this model and the Coulomb gas in the Kostelitz-Thoulcss phase 
(/3 > Stt), go along with the non-existence of any proof for the analyticity of the 
pressure in the Coulomb gas. Indirect arguments are attempted in [62], [63] and 
[64] . We believe that after a suitable regularization of the Coulomb potential at 
short distances to assure stability, we can fit the problem into the framework of 
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the model discussed above and get an estimate of the coefficients of the Mayer 
series through our formula for the nth derivative of the pressure. 
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